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Abstract 

BES data on the a pole are refitted taking into account new information on coupling of a 
to KK and 7777. The fit also includes Cern-Munich data on tttt elastic phases shifts and K e ^ 
data, and gives a pole position of 500 ± 30 — i(264 ± 30) MeV. There is a clear discrepancy 
with the a pole position recently predicted by Caprini et al. using the Roy equation. This 
discrepancy may be explained naturally by uncertainties arising from inelasticity in KK 
and r/77 channels and mixing between a and /o(980). Adding freedom to accomodate these 
uncertainties gives an optimum compromise with a pole position of 472 ± 30 — £(271 ± 30) 
MeV. 



PACS: 13.25.-k, 13.75.Lb, 14.40.Cs, 14.40.Ev. 
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\Q BES data on J/$ -> um+ir display a conspicuous low mass tttt peak due to the a pole [1]. It 
was observed less clearly in earlier DM2 [2] and E791 data [3]. The BES data are reproduced in 
Fig. 1. The band along the upper right-hand edge of the Dalitz plot, Fig. 1(a), is due to the a 
1 pole. There is a clear peak in the nn mass projection of Fig. 1(b) at ~ 500 MeV; the fitted a 
^■contribution is shown by the full histogram of Fig. 1(d). Other large contributions to the data 
^ arise from / 2 (1270) and 61 (1235), which appears in the um mass projection of Fig. 1(c). 

! At the time of the BES analysis, little was known about the coupling of a to KK and r/r/, and 
^ jthese channels were omitted from the fit to the data. Since then, couplings to A" A" and r/r/ have 
S-i been determined by fitting (a) all available data on mi — > KK and r/r/, (b) Kloe and Novosibirsk 
data on <fi — > 7(7r°7r°) [4]. All these data agree on a substantial coupling of o to AA. The first 
objective of the present work is to refit the J/^ — > uj-k-k data including this coupling, following 
exactly the procedure of the BES publication. The outcome is to move the a pole position from 
541 ± 39 - i(252 ± 41) MeV reported in the BES paper to 500 ± 30 - i(264 ± 30) MeV. 

Meanwhile, Caprini, Colangelo and Leutwyler (denoted hereafter as CCL) have made a pre- 
diction of the a pole position [5]. Their calculation is based on the Roy equation for irn elastic 
scattering [6], which embodies constraints of analyticity, unitarity and crossing symmetry. This 
approach has the merit of including driving forces from the left-hand cut due to exchange of p, f% 
and a. They also apply tight constraints from Chiral Perturbation Theory (ChPT) on the S-wave 
scattering lengths a and for isospins and 2. Experiment alone determines only the upper 
side of the pole well; this has led to speculation that fits without a pole might succeed in fitting 
the data [7]. CCL's use of crossing symmetry and analyticity gives a precise determination of the 
magnitude and phase of the S-wave amplitude on the lower side of the pole and leaves no possible 
doubt about its existence. They clarify the fact that tttt dynamics are fundamental to creating 
the pole. They quote rather small errors: M a = 441 tg 6 — 2(272^12.5) MeV. There is a rather large 
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Figure 1: BES data of Ref. [1]. (a) Dalitz plot; (b) tttt mass projection: the upper histogram shows 
fit (ii) of Table 1, the lower histogram shows experimental background; (c) ujtt mass projection; 
full histograms are as in (b) and the dashed histogram shows the coherent sum of both bi(1235)n 
contributions; (d) tot mass projections from a (full curve) and spin 2 (dashed). 
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discrepancy between this prediction and BES data. A second objective of the present work is to 
trace the origin of this discrepancy. In outline, what emerges is as follows. 

From CCL results for T 11 (s), the T-matrix for irn elastic scattering, one can predict what 
should appear in BES data. The prediction will be shown below on Fig. 4 by the chain curve. 
From a glance at this figure, one sees a significant disagreement with the experimental points, 
which are deduced directly from BES data. The questions which arise are as follows: 

• (i) are the hypotheses used in the BES analysis wrong or questionable? 

• (ii) is something missing from the calculation with the Roy equation? 

• (iii) can the calculation of CCL be fine-tuned in order to come into agreement with the 
experimental data? 

The discrepancy lies in the mass range 550 to 950 MeV. In this range, the analysis is sensitive 
to assumptions about inelasticities due to KK and 7777 channels. These are not known with 
sufficient accuracy at present and allow freedom in the analytic continuation of couplings to these 
channels below their thresholds. The amplitude for KK — > tttt goes to zero at the KK threshold; 
it also has an Adler zero at s ~ 0.5m^. In between, it has a peak near 500 MeV providing a 
natural explanation of the additional peaking required by BES data. This can readily explain the 
discrepancy, as shown by the full curve on Fig. 4. It may be regarded as a fine-tuning of the 
solution of the Roy equation. It will be shown that changes to nn phase shifts up to 750 MeV are 
< 1.2°, well below experimental errors. In summary, the Roy equation and Chiral Perturbation 
Theory provide the best description of tttt scattering near threshold (and below), while the BES 
data provide the best view of the upper side of the a pole. 

The layout of the paper is as follows. Section 2 introduces the equations. Section 3 discusses 
the prediction of CCL and subsection 3.1 explains what is missing from their Roy solution. Section 
4 describes fits to BES data. Section 5 discusses possible alternative explanations of BES data 
and Section 6 summarises conclusions. 

2 Equations 

The 7in elastic amplitude may be written 

T n (s) = ( V e™-l)/m = N(s)/D(s), (1) 

where N(s) is real and describes the left-hand cut; D(s) describes the right-hand cut. The nu- 
merator contains an Adler zero at s = sa — 0.41m^.: 

N(s) = (s-s A )f(s), (2) 

where f(s) varies slowly with s. The standard relation to the S-matrix is Tn = (Sn — l)/2i; 
Tu has contributions from / (980), / (1370) and / (1500) as well as a. For elastic scattering, 
these contributions are combined by multiplying their S-matrices, to satisfy unitarity; this implies 
adding phases rather than amplitudes. 

In the BES publication, it was shown that their data may be fitted with the same 1/D(s) 
for the a component as K e4 data [8] and Cern-Munich data [9] within errors. The phase of the 
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nn S-wave in BES data has the same variation with s as the a component in elastic scattering 
within experimental errors of ~ 3.5° from 450 to 950 MeV [10]. The BES data may be fitted by 
an amplitude A/D(s); here, A is taken as a constant, since the left-hand cut of J/ty — > uoa is very 
distant. Possible doubts about this assumption will be discussed later. Channels nn, KK, rjr] and 
An will be labelled 1 to 4. The parametrisation of the a is given by 
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D(s) is the denominator of eqn. (3). Below the KK threshold, Im D = —N(s) and tan S = 
—Im D/Re D, from which one can deduce D(s). Since D(s) is analytic, it is subject to a 
normalisation uncertainty which is constant within the errors of T n (s), in practice 2-3%. 

The function ji(s) is obtained from a dispersion integral over the phase space factor p 7r7r (s) = 
An important point about z(s) 



1 — Am\js. An important point about z(s) is that it is well behaved at s = and eliminates 
the singularity in p(s). It makes the treatment of the -kit channel fully analytic - an improvement 
on earlier work. In the An channel, a dispersion integral is in principle required but is small below 
1 GeV and can be absorbed into the fit to g 2 {s). In eqn. (11), An phase space is approximated 
empirically by a combination of pp and aa phase space [11] (with s in GeV 2 ); T(47r) is set to zero 
for s < 16m 2 . The effect of the An channel on the a pole is only 2 MeV and is not an issue. 

It is assumed that the Adler zero is a feature of the full nn amplitude. The factor (s — 
s A )/(M 2 — sa) of eqn. (4) introduces this Adler zero explicitly. Eqn. (5) is an empirical form 
used earlier in fitting BES data; b±, b 2 and A are fitted constants. The factors s/M 2 in T 2 and 
Ts of eqns. (8) and (9) approximate the Adler zeros closely at s = 0.5m 2 , 0.5m 2 K and 0.5m 2 and 
remove the square root singularity in p 2 and p 3 of eqns. (8) and (9). 
g 2 KK and g 2 have been fitted to data on nn — > KK and rm and on < 
also determine olrk = OL m = 1.3 GeV 2 above the thresholds. 

A detail is that the factor (s- s A ) / (M 2 - s A ) is also used for r\, T 2 and T 3 of / (980), / (1370) 
and /o(1500). Otherwise, parameters of /o(980) are taken from BES data on J — > <pn + n~ and 
<f)K + K~ [12]. The /o(1500) and /o(1370) are fitted with Flatte formulae and parameters given in 
Ref. [4]. The combined contribution of / (980), / (1370) and / (1500) to the scattering length is 
8%. 

An important question is how to parametrise the continuation of KK — > nn and r\r\ — > nn 
amplitudes below their thresholds, i.e. in K- matrix notation the elements K\ 2 and K\o,. They can 
in principle be determined from dispersion relations for each amplitude. The factors p 2 and ps 



The factors 0.6 and 0.2 for 
) -> 7(7r°7r°) [4]. These fits 
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in eqns. (8) and (9) are kinematic factors. Below the thresholds, they are continued analytically 

as i\J<imf — s. If p 2 is factored out of the KK — > nir amplitude, there is a step at threshold 
in the imaginary part of the surviving amplitude. This step leads to a cusp in the real part of 
the amplitude at threshold, i.e. a change of slope in eqns. (8) and (9). An evaluation of the 
dispersion integral generates a result for the real part of the amplitude below threshold close to 
an exponential falling as exp[— a (4m? — s)]. 




Figure 2: (a) the magnitude of the KK — > tttt amplitude from Biittiker et al. [13]; the result from 
my fit to Kloe data (full curve) and the a contribution alone (dashed); (c) fit to Kloe — > 7(7r 7r°) 
data after background subtraction; (d) contributions from / (980) (full curve), a (dashed) and 
interference (dotted). 



A calculation of Ki 2 has been made along these lines by Biittiker, Descotes-Genon and Mous- 
sallam using the Roy equations for nK — > nK and ixix — > KK [13]. Their result is reproduced in 
Fig. 2(a). The peak is due to / (980) and the low mass tail comes from the analytic continuation 
of g — > KK. Using their K\ 2 in fitting BES data does lead to effects with the right trend, but 
not with sufficient accuracy to make a good prediction. They remark that their calculation is 
uncertain because of discrepancies between available sets of data on ttti — > KK. My estimate is 
that errors of Caprini et al from the dispersion relation are at least a factor 2 too small. 

My conclusion is that K 12 and K^ presently need to be fitted empirically. There is however 
experimental information which helps decide an appropriate parametrisation: data on — > 7(7r°7r°) 
from Kloe [14]. In Ref. [4], those data are fitted using the standard KK loop model of Achasov 
and Ivanchenko [15]. Parameters of /o(980) are fixed to those determined by BES data [12]. The 
fit requires a substantial additional amplitude for KK — > nn through the a. Empirically this 
contribution is well fitted with the exponential of eqn. (8) with a = 2.1 GeV -2 below the KK 
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threshold and an Adler zero at s = 0.5m 2 K . [The same a is assumed for r\r\ — > -kit.] Fig. 2(c) 
reproduces from Ref. [4] the fit to Kloe data; Fig. 2(d) shows the / (980) and a components and 
the interference between them. Fig. 2(b) compares my fit with that of Biittiker et al.; the dashed 
curve shows the a contribution. 

3 The Roy solution of CCL 

CCL make a prediction of the tttt — > tttt S-wave amplitude and deduce the a pole position from 
it. The inputs to their calculation are [5,16,17]: 

i) the Roy equation, which accounts for both left and right-hand cuts in tttt elastic scattering; 

ii) the precise lineshape of p(770) from e + e~ — > tt + tt~ data; 

iii) predicted value from ChPT for S-wave scattering lengths a = (0.220 ± 0.005)m~ 1 , a 2 = 
(-0.0444 iO.OOlO)™; 1 [17]; 

iv) the tttt phase shift at 800 MeV with an error of i^ deg; 

v) my elasticity parameters 77 above the KK threshold; 

vi) minor dispersive corrections for masses above 1.15 GeV, their matching point. 
Note that experimental phase shifts are not fitted except for constraint (iv). 

CCL have kindly supplied a tabulation of their Tn(s). My first step is to reproduce this solution 
using eqns. (1)— (11). This is simply a fit to their fit, i.e. an explicit algebraic parametrisation. 

They do not explicitly separate a and / (980). This leads to uncertainty in how the / (980) 
is being fitted. In the vicinity of the a pole, an issue is the magnitude of the low mass tail of 
/o(980), which affects what is left as the remaining a amplitude. CCL omit / (1370) and /o(1500) 
contributions and a — > 4tt. 

The fit shown in column (i) of Table 1 is made in two steps. Firstly a fit to their / (980) is made 
from 800 to 1150 MeV. It requires M = 0.970 GeV, g\ = 0.146 GeV 2 , somewhat smaller than 
0.165 GeV 2 from BES data on J/* -> 0/ o (98O) [12]. This step also reveals that their a -> KK 
and rji] components fall from these thresholds at least as fast as s 2 25 ; their contributions below 
800 MeV are very small. In the second step, the mass range below 800 MeV is refitted, in order 
to minimise the sensitivity of a parameters to the f mass region. Empirically, this second step 
can reproduce CCL phases only if (a) a — > KK and r\r\ contributions are omitted and (b) the 
fo — > tttt amplitude is multiplied by a factor (s/M 2 ) n , where n optimises at 1.55. Their phases 
are then reproduced everywhere up to 750 MeV within errors of 0.12°, as shown by the full curve 
of Fig. 3. This figure shows fitted values minus the phase shifts of CCL for the full tttt elastic 
amplitude. The pole position is 449-i269 MeV, 8 MeV higher in mass than the CCL pole. 

The essential conclusion is that their contributions from / (980) — > tttt, KK — > tttt and rjrj — > 
tttt are cut off very sharply at high mass, and the whole of the tttt amplitude below about 500 
MeV is being fitted by the a alone. This is hardly surprising. It is well known that /o(980) is 
driven largely by forces in the KK channel, not by the left-hand cut of the tttt channel. This 
question will be discussed more fully below. For the moment, it is sufficient to remark that the 
pole immediately moves up to 467 MeV if the / (980) — > tttt amplitude is multiplied only by the 
factor s/M 2 of eqns. (8) and (9) and the KK — > tttt and rj — > tttt amplitudes of the a are treated 
with the falling exponential fitted to Kloe data. 

Fig. 4 displays values of |l/-D(s)| 2 from BES data as points with errors, normalised to 1 at 
their peak; they are obtained by dividing out 3-body phase space from data of Fig. 1(b). Note 
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(i) 


(ii) 


(iii) 


M (GeV) 


1.038 


0.958 


0.953 


61 (GeV) 


1.082 


1.201 


1.302 


b 2 (GeV- 1 ) 
A (GeV 2 ) 


-0.016 


0.684 


0.340 


1.179 


2.803 


2.426 


S-47T (GeV) 





0.014 


0.011 


Pole (GeV) 


0.449 - i0.271 


0.500 - i0.264 


0.472 -i0.271 



Table 1: Parameters of fits discussed in the text. 
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Figure 3: Fits (i)-(iii) of Table 1 minus the predicted phases of CCL. The full curve is fit (i) 
to CCL phases, the dashed curve is fit (iii), the compromise fit between CCL phases and BES 
data, the chain curve is fit (ii) to K e4 and Cern-Munich phases and BES data; points with errors 
illustrate departures from CCL phases of K e4 , Cern-Munich data and a point at the K mass. 




that the normalisation is arbitrary - it is the s-dependence which matters. There is a significant 
disagreement with the prediction from fit (i) to CCL phase shifts (chain curve); this discrepancy 
is far larger than the 2-3% errors in deducing D(s). The essential point is that BES data fall more 
rapidly from 600 to 950 MeV than CCL's result. Fig. 5 shows the poor fit to BES data using their 
D(s). Note that the discrepancy is not a question of the extrapolation of amplitudes to the pole. 
There is a direct conflict for physical values of s between CCL and the fit to BES data. However, 
bearing in mind that this is a theoretical prediction, the closeness to data is still remarkable. 

3.1 Missing elements in the Roy solution 

The Roy equation is in principle exact. How is it possible to modify the solution? In Nature, the 
Roy equation 'knows' about coupling of ixix to KK and 7777, as well as the /o(980), /o(1370) and 
/o(1500) resonances. If tttt phase shifts were known with errors of a small fraction of a degree, 
it would be possible to deduce these resonances and their couplings. However, that is not the 
present situation. In reality, there are subtle features in the processes irn — > KK and r\r\ arising 
from meson exchanges and also from mixing between a and /q. 
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M (GeV) 

Figure 4: 1 1 / -D(s) | 2 from BES data ; the chain curve is fit (i) to T\\ values of CCL. The dashed 
and full curves are from fits (ii) and (iii) of Table 1. 




Figure 5: Fit to BES data using 1/D(s) derived from CCL phases. 



Although the a pole appears remote from the KK threshold, one must remember that the 
phase shift it produces reaches 90° close to the position of /o(980), so multiple scattering is a 
maximum there. This includes terms of the form a — > nn (or KK) — > /o(980) leading to mixing. 
Anisovich, Anisovich and Sarantsev show [18] that this mixing obeys the Breit-Rabi equation, so 
a and /o(980) behave as a pair of coupled oscillators. 

A recent paper with van Beveren, Rupp and Kleefeld [19] shows that the nonet of a, k, ao(980) 
and /o(980) may be generated by coupling of mesons to a qq loop. The a is generated by the nn 
loop and the /o(980) by the ss loop. Without mixing between a and /o, one can only account for 
< 10% of the observed magnitude of / (980) — > nn. Substantial mixing is required to produce 
the observed nn width. Using their programme, I have tried varying the mixing angle and its 
s-dependence to examine perturbations of the nn amplitude. Since the Schrodinger equation 
is solved, the amplitudes are fully analytic. In the mass range 700-1150 MeV, one sees subtle 
correlated changes in nn phases and the inelasticity parameter rj. These are not present in the 
calculation of CCL. The Roy equation is simply a dispersion relation between real and imaginary 
parts of the nn — > nn amplitude. Its solution for the real part is only as good as the input for 
the imaginary part. They derive results for / (980), KK — > nn and nn — > nn from my inelasticity 
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parameters of Ref. [4]; but those are based upon simple Flatte parametrisations of f and a, 
without any dynamics due to mesons exchanges and mixing. 

Calculations in Ref. [19] show that the f pole is anchored to the KK threshold. What 
happens as the mixing is varied is that an f — > tttt amplitude reaching down to the tttt threshold 
is generated. It is questionable exactly what its precise mass dependence will be. The best 
conjecture one can make at present is that it will contain the Adler zero of the full tttt amplitude, 
but otherwise behaves like a normal resonance. That is the conjecture adopted here. Calculations 
with the model of van Beveren et al show that the a pole can be affected by the mixing with 
/o(980) by up to 50 MeV. The reason for this is that the mixing alters the s-dependence of the tttt 
amplitude. The pole position is determined by the Cauchy-Riemann relations as one moves off the 
real s-axis. It lies far below the mass at which the tttt phase shift reaches 90°. Because of this long 
lever arm, even a small change of the slope of tttt phases v. mass can move the pole a surprisingly 
long way. This is why including an / (980) contribution making up 7% of the scattering length 
moves the pole position up to 467 MeV. 

This degree of freedom is missing from the calculation of CCL. My conclusion is that it is 
legitimate to introduce into the fit small systematic perturbations in eqns. (4), (8) and (9) arising 
from inelasticities in KK and rjr] channels and also mixing between a and fa. This secures 
agreement with BES data with only small changes in tttt phase shifts. 

4 Refit of BES data 

In the next fit, BES data are refitted together with K e 4 and Cern-Munich data as in the BES pub- 
lication, but using the new equations of Section 2. The small /o(1370) and /o(1500) contributions 
are included in this fit (and the next). The scattering length a and the effective range are fitted 
with the errors quoted by CCL. Parameters are shown in column (ii) of Table 1. 

Resulting errors decrease compared with the BES publication because including a — > KK and 
rjr) produces cusps at these thresholds and removes the requirement for the /o(980) contribution 
completely; the errors of ±30 MeV in real and imaginary parts of the pole position are taken from 
a range of fits to BES data using a variety of fitting functions going beyond those used here; they 
include systematic errors in data. The shift in the pole arises partly from the use of ji(s) of eqn. 
(6) but mostly from including KK and 7777 channels. There is a distinct improvement in the fit to 
BES data. This fit is used in making Fig. 1. It is a valuable feature of the BES data that the 
/o(980) contribution is negligible, giving an unimpeded view of the a pole. 

The fit to phases is shown by the chain curve of Fig. 3 and the fit to BES data by the dashed 
curve of Fig. 4. The fitted value of a falls to O.lSdm^ 1 , i.e. Q.2a below the ChPT prediction. 
The question is whether to believe the ChPT prediction or the experimentally fitted scattering 
length. The ChPT prediction was made using the Roy equation which gives the low mass of 
CCL's a pole. It is obvious that there will be a correlation between ao and the pole position: as 
the pole moves away from the tttt threshold, the scattering length will naturally decrease. If the 
scattering length is forced upwards, one finds an almost linear relation between the pole mass and 
the scattering length. To reach a value ao = O^lSm^ 1 requires a pole mass of 472 MeV. The \ 2 
of the fit increases substantially, mostly because of the K e 4 and Cern-Munich data; in particular, 
the 382 MeV point of K e4: data pulls the scattering length down. The BES data at low mass also 
favour a small scattering length, but are in a mass range where the available phase space is cutting 
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off the signal. 

It must be said that this fit is not using information from the Roy equation. It is made 
empirically to data above threshold. In the next step, the information from the Roy equation is 
incorporated as closely as possible. 

4.1 A compromise solution 

In order to constrain the fit to data as closely as possible to the Roy equation, a combined fit is 
made to the T\\ of CCL up to 1.15 GeV, and data from BES, and Cern- Munich. Errors are 
assigned to CCL phases rising linearly with lab kinetic energy from zero at threshold to 3° at 750 
MeV. This gives maximum weight to results from the Roy equation near threshold, but allows 
flexibility in the effects of / (980), and continuations of effects of KK and rjrj inelasticities into 
the mass range 600 to 1000 MeV. The pole moves to (472 - i 271) MeV and parameters are given 
in column (iii) of Table 1. The scattering length is 0.215m" 1 . The fit to BES data is shown by 
the full curve of Fig. 4; it is marginally poorer in terms of total \ 2 but is obviously acceptable. 
However, the fit resists strongly any attempt to move the mass of the pole any lower and the 
scattering length correspondingly higher. 

The fit to phases is shown by the dashed curve of Fig. 3. Up to 750 MeV, the difference 
between this fit and CCL phases is only 1.2° at 470 MeV. This illustrates how difficult it is to 
deduce the a pole position from elastic data. The difference of 1.65° at 865 MeV is well within 
errors of data. 

Strictly speaking, the small perturbation to the a amplitude on the right-hand cut induces, via 
crossing, a small perturbation to the left-hand cut. However, the a contribution on the left-hand 
cut has a small isospin coefficient, and I have checked that in practice the perturbation is smaller 
than errors in the large contributions from p and j'2 exchange. 

What scattering length should be adopted? The lowest order prediction from Weinberg's 
current algebra [20] is 0.16m" 1 . In second order ChPT it rises to 0.20m,,: 1 and then 0.215m" 1 
at fourth order. The prediction of (0.220 ± 0.005)m n T 1 at sixth order by Colangelo et al. [17] 
takes account of the unitarity branch cut at threshold; however, it does depend on using the 
Roy equations. It seems safe to constrain the scattering length to be at least 0.215m" 1 , so the 
conclusion is that fit (iii) is currently the best compromise. The change to the pole mass predicted 
by CCL is modest, but twice the error they quote. 

The full and dashed curves of Fig. 4 are very similar. This is because they differ only in the 
scattering length ao- An interesting result is that, for all three fits, the effective range does not 
change significantly from the CCL value. In lowest order ChPT, the effective range is proportional 
to 

2a - 5a 2 = 3ml/ An F^, ( 12 ) 

as shown by Weinberg [19], so this relation is accurately consistent with all three fits. In this 
relation, F n is the pion decay constant. 

5 Possible ambiguities in fitting BES data 

Although the effect of systematic changes in inelasticity over the mass range from threshold to 1.15 
GeV provides a natural resolution of the discrepancy between BES data and the CCL prediction, 
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alternatives have been suggested. Two points have appeared in recent preprints. Firstly, Wu and 
Zou [21] remark that the width fitted to 61 (1235) is 195 MeV compared with the value 142 ±9 MeV 
of the Particle Data Group [22]. This discrepancy was also reported by DM2 [2]. Wu and Zou 
suggest that the strong process J/^J — > pn, followed by p — > ujtt may contribute. In the original 
BES work, the possibility p(770) — > ojix was tried and gave little improvement and, at maximum, 
a contribution of 2% (intensity) of the data. In any case, the effect lies close to the vertical and 
horizontal b\ bands of Fig. 1(a), particularly the region where they cross near the bottom left-hand 
corner of the Dalitz plot. This is remote from the a band; changes to interferences between a and 
61, p or p' have negligible effect on a parameters. 

Secondly, Caprini [23] suggests that triangle graphs due to b\ — > uj-k, followed by nn rescattering 
will introduce effects beyond the isobar model. Although this is true, it is known that such effects 
vary logarithmically over the Dalitz plot. There is no obvious reason why they should introduce 
a rapidly varying effect close to the right-hand edge of the Dalitz plot. 

Thirdly, could there be some form of background in J/ty — > uoal If such a background is 
included as a quadratic function of s, the discrepancy with CCL persists. The reason is that 
BES data determine the nn S-wave amplitude accurately above the / 2 (1270) as well as below it 
and limit it to small values; this severely limits the form of any background to something close 
to the observed peak. In particular, a conventional form factor exp(— k 2 R 2 /6) at the vertex for 
J/ty — > uj<7 (where k is u momentum in the J/\& rest frame) gives no improvement, since it requires 
negative R 2 . So a background does not provide a simple escape route. 

A fourth possibility has been raised in discussions with CCL. This is that the form factor has a 
zero somewhere above 1 GeV. Extrapolating the full curve of Fig. 4, there might in principle be a 
zero in the mass range around 1.3 GeV. This would be obscured in the BES data on J/$ ^ ujtxtx 
by the /2(1270). However, there are also data on J/$ ^ uK + K~ where the a — > KK amplitude 
is clearly required [24]. The / 2 (1270) — > KK contribution in those data is sufficiently small that a 
zero in the a — > KK amplitude can be ruled out definitively up to ~ 1.6 GeV. At that mass, the 
required radius of interaction would be unreasonably large, > 0.8 fm, and the a pole region would 
be seriously distorted by the form factor. Data on production of the k pole in <//*& — > K*(890)K7r 
require an RMS radius < 0.38 fm with 95% confidence [25]. 

6 Discussion and conclusions 

There is a significant conflict in Fig. 4 between the a pole of CCL and BES data. In my view, the 
discrepancy needs explanation and the BES data should be taken at face value. The prescription 
adopted here in eqns. (4), (8) and (9) gives a natural improvement in fitting BES data without 
disturbing elastic phases up to 750 MeV by more than 1.2° in fit (iii). This illustrates the ease 
with which the prediction of CCL may be modified to fit BES data. 

The strength of the BES data is that a peak is clearly visible. There is no significant /o(980) 
signal. The BES data are therefore free of uncertainties about mixing between a and / (980). The 
weakness of elastic scattering data is that there is no visible peak which can be checked, and there 
is a significant / (980) amplitude which must be separated. The BES data provide a better view 
of the upper side of the a pole than elastic data, where the low mass tail of /o(980) is uncertain. 
The CCL calculation provides a better view of its lower side, where constraints from ChPT are 
valuable. 



11 



A more ambitious approach, beyond the scope of present work, would be a solution of coupled 
channel Roy equations for nn, KK and rjrj channels, including the dynamics driving / (980) and 
a — > KK and rji]. However, from experience with the model of van Beveren et al [19], it is likely 
that uncertainties are presently too large to give a definitive prediction of the dynamics of /o(980) 
and its delicate mixing with the a. 

I am grateful to Prof. H. Leutwyler for extensive discussions of the details of the calculation 
of Ref. [1]. I am also indebted to Profs. G. Rupp and E. van Beveren for use of their programme 
and Dr. F. Kleefeld for illuminating discussions. 
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